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We examine the effect of a nonmagnetic impurity in a two-dimensional spin liquid in the spin- 
gap phase, employing a drone-fermion representation of spin-1/2 operators. The properties of the 
local moment induced in the vicinity of the impurity are investigated and an expression for the 
nuclear-magnetic-resonance Knight shift is derived, which we compare with experimental results. 
Introducing a second impurity into the spin liquid an antiferromagnetic interaction between the 
moments is found when the two impurities are located on different sublattices. The presence of 
many impurities leads to a screening of this interaction as is shown by means of a coherent-potential 
approximation. Further, the Kondo screening of an impurity-induced local spin by charge carriers 
is discussed. 
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I. INTRODUCTION 

Substitution of Cu ions in the conduction planes of 
high-Tc cuprates by different nonmagnetic ions presents 
an important experimental tool for the study of the 
metallic state of these strongly correlated systems. Un- 
usual effects have been revealed especially when these 
materials were doped with Zn, Al, or Ga. Among 
these interesting features is the appearance of a mag- 
netic moment derived from the observation oLa Curie- 
Weiss behavior of the magnetic susceptibility.El Further 
studies have been performed by electron paramagnetic 
resonancea (EPR) and nuclear magnetic resonancetTEl 
(NMR) experiments. Measurements of the NMR 
Knight shift indicate that the impurity-induced local mo- 
ments i"fij4de predominantly on Cu sites neighboring the 
dopant .Bi3 It has been argued that the appearance of 
these locaLmoments can account for a significant reduc- 
tion of Tc.B 

On the theoretical side, impurity-induced moments 
have been studied for a variety of quantum disor- 
dered spin systems, includins, r-jfor instance, spin lad- 
ders and|-S|pin-Peierls systemsJ3'u as well as underdoped 
cuprates Jj'E2l A common feature of these systems is the 
existence of a (pseudo-) gap in the spin-excitation spec- 
trum. This spin gap is related to the singlet-pairing cor- 
relations in the ground state, which can be described in 
terms of resonance vaWe bondsEl (RVB) or a valence- 
bond soHd (VBS) state.EHl Within this picture the appear- 
ance of a Curie-type susceptibility can be understood as 
an unpaired spin that is generated by the substitution of 



a spin site by a vacancy. 

In the present paper we investigate the effect of 
nonmagnetic impurities on the local magnetic proper- 
ties of weakly doped cuprates. In order to describe 
the spin degrees of freedom of these systems we start 
with the Heisenberg antiferromagnet on a square lat- 
tice, employing_.a drone-fermion representation of spin- 
1/2 operators.!!^ A mean-field decoupling of the Hamil- 
tonian leads to solutions corresponding to a RVB state 
in the fiux phaseoEJ that exhibits a pseudogap in the 
spin-excitation spectrum. A single spin vacancy intro- 
duced into the spin liquid produces a local moment which 
is predominantly located on the sites close to the va- 
cancy. We derive an expression for the impurity-induced 
NMR Knight shift and fit it to experimental data. In- 
troducing a second vacancy into the spin liquid leads to 
an antiferromagnetic interaction between the two mo- 
ments. In the presence of many impurities this interac- 
tion is screened as is shown by using a coherent-potential 
approximation (CPA). We further analyze the effect of 
nonmagnetic impurities on the relaxation of the nuclear 
magnetic moment in terms of local spin fluctuations in- 
duced by the vacancies. Finally, we discuss qualitatively 
the low-temperature Kondo screening of the impurity- 
induced spins by charge carriers. 

II. THE MODEL 

We start from the Hamiltonian of a two-dimensional 
square lattice spin-1/2 Heisenberg antiferromagnet 
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H 



(1) 



where ( ) indicates summation over pairs of nearest- 
neighbor sites. The spin operators in the Hamiltonian 
(|l|) are commonly expressed in terms of pseudo-fermions 
which requires us to impose the restriction that each site 
is occupied by exactly one pseudoparticle. Since it is 
rather difficult to account for these local constraints one 
usually treats them onlyjiin the average. Here we will fol- 
low a different approachtS by employing a dipne-fermion 
description for spin- 1/2 operators given byll3 



fiXi^ — Xifii — fifi ■ 



(2) 



where Sf" — {Sf ± iSf)/y/2. The operator fj creates a 
spinless fermion at site i; the presence (absence) of an /- 
fermion corresponds to an up (down) spin state. The real 
drone-fermion operator Xi = xl with commutation rules 



[Xr^Xj] 



Sij is needed to provide the proper commuta- 



tion rules for spin operators on different sites. Expressing 
the Hamiltonian (|^) in terms of the representation (^) one 
obtains 



H = J 



E{(//xa,/,+H.c, 

iv) 



(3) 



Mean-field solutions to the Hamiltonian (Q) corre- 
sponding to the RVB states of conventional fermion 
mean-field theory are found by employing a Hartree-Fock 
factorization with bond parameters iA.ij ~ (xjXi) ^^'^ 
iAij ~ {fjfi) defined for nearest-neighbor pairs, yielding 
the following mean-field Hamiltonian:0 



i?MF — J 



iv) 



- ^A,JX^XJ + H.c. ) + 2A.yAy- + A^J. 



(4) 



For spin rotational symmetry of the ground state A^ — 
Aij holds. In the following we assume the order param- 
eter to depend only on relative site indices: A^ = A^ 
with 6 £ {x, —x,y, —y}, where x and y are the lattice 
unit vectors. 

Dividing the square lattice into two sublattices A and 
B, we define i £ A and j E B. The mean-field Hamilto- 
nian (^) can then be diagonalized by expressing the / and 
X operators in momentum representation and employing 
the canonical transformation 



= (flfc + bk) /V2, 
fk = [o-k - bk) ex.p{iipk)/V2, 
xt = {ck + dk)/V2, 
Xk = (ck - dk) exp{iLpk)/V2, 



(5) 
(6) 



where exp(i(pfc) = 
^Z^a '^s exp{-ik5) 



Ak\ with the definition Ak = 
Thus, the Hamiltonian (H) becomes 



-ffMF = S.k (-alflfc + blbk 

fceMBZ 

+ ^ Cfe (^-clck + dldk 

fceMBZ/2 



(7) 



with ^k = 2J|Afe|. The k summations extend over the 
full and half magnetic Brillouin zone, respectively, and 
the constant terms in Eq. (^ have been omitted. 

The fermionic spectrum E{k) — ±^fc is determined by 
the symmetry exhibited by the bond parameter A^. In 
a phase of mixed symmetry defined by 



A Ai = -A_£ = Ay = A^y 



one finds 



Cfc = \l sin"^ kx + cos^ ky , 



(8) 



(9) 



where D = 4-\/2JA is the half-width of the band and 
the lattice constant has been set equal to unity. At low 
energies the density of states is determined by the neigh- 
borhood of the two isolated roots of Eq. (p|) and exhibits 
a V-shaped pseudogap p°(w) = 2|aj|/(7rlP) correspond- 
ing to the one found in th e | flii x phase of conventional 
fermion mean- field theories .Ejlla 

The bandwidth parameter D in Eq. (|^) is obtained 
from the self-consistency relation 



N ^ \2V2tJ 



(10) 



where rjk — (siia^ k^ + cos^fcy)^/^, N is the number of 
sites, and fc^ = 1- In the limit of T ^ Eq. ( |lO| ) reduces 
to 



Dn 



N 



E 



Vk, 



(11) 



yielding a numerical value of Dq 1.355 J. The cor- 
rection SD{T) to Do for small but finite temperatures is 
given by 



dD{T) 



32J / T 



(12) 



In the following we study the effect of a spin va- 
cancy located at site (sublattice A) on the local mag- 
netic properties of the system described by Hamiltonian 
(0). We simulate this spin defect by decouphng site 
from the rest of the system. This is done in two steps: 
First, the drone-fermion bond parameter connecting the 
impurity site with its nearest neighbors is set to zero, 
iAo^s — ixsXo) = 0. This decouples the drone fermions 
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of site from the rest of the system. Then, a local chem- 
ical potential Xq oo acting on site is added to Hamil- 
tonian 



(13) 



It induces an empty site at i? = with respect to 
the spinon / operators and suppresses the low-energy 
/-fermion degrees of freedom at this site. We will show 
below that the presence of the local potential also leads to 
a vanishing spinon bond parameter connecting the impu- 
rity site with its nearest neighbors, iAo^s = {fsfo) — 0, 
and therefore completely decouples the impurity. Fur- 
ther, one infers Ah,h+5 — 'Ar,r+s for all R £ A 
which manifests the equivalent treatment of / and drone 
fermions as required by spin rotational symmetry. 



III. LOCAL MOMENT 

In this section we study the uniform properties and 
the spatial distribution of a local moment induced by a 
spin vacancy. We first calculate the impurity contribu- 
tions to the spinon density of states and the uniform spin 
susceptibility by means of the uniform Green's function 
G(*c.„) = -l/Nj:^{Trfr{T)fm)^„ given by 



(14) 



with fermionic frequencies w„ = (2n + l)7rT'. The Green's 
function in the absence of the impurity G^{iujn) and the 
correction SG{iLUn) due to the impurity are given by 

G''(*c^„) = ^^5°(*'^n;fe), 

k,i^ 

SG{iujn) = X! 5°(«^n; k)T{iuj,i)gl{iuJn; k) 



(15) 



with g°(zw„;fc) = l/{tuj.n - (-1)"^), 1^ G {1,2}. The T 
matrix in Eq. (15) describes successive scattering at the 
vacancy and is obtained from usual scattering theory; in 
the limit Aq — > oo it is T(iw„) — — 1/G°(ia;„), yielding 



d 



(5G(^c^„) = -— lnG°(zc^„) 
mujn 



Evaluating Eq. (|l^) one obtains 



i5G(iw„) 



1 - In" 



D 



(16) 



(17) 



From Eq. (|T^) we calculate the impurity correction to 
the density of states Sp{lu) ~ — l/7rIm(5G((jj + zO+): 



2|c^| {^/2f +\n^D/\u;\)' 



(18) 



The first and second terms on the right-hand side of 
Eq. (nsh describe, respectively, the formation of a spinon 



bound state at the Fermi surface and the destruction of 
a spinon singlet in the spin background. Integrating Eq. 
( |T8| ) over all occupied states one finds a reduction in the 
number of particles by 1/2. This is in agreement with the 
fact that the local chemical potential Aq induces an empty 
site at i? = with respect to the / fermions, while in the 
absence of Ao the average /-fermion occupation number 
is 1/2. 

With the help of Eq. (|l^) we evaluate the impurity 
contribution 5x{T) to the uniform spin susceptibility by 
employing the relation 



5x{T)^TY,J-m^^n). 

^ — ' OlUJn 



For T J one obtains 
1 



5x{T) 



4r V T 



(19) 



(20) 



The formation of a spinon bound state is thus found to 
lead to a Curie-like spin susceptibility of a free spin 1/2 
with logarithmic correcticp. We note that conventional 
fermion mean-field theoryll3 only reproduces one half of 
the above result. This shortcoming has its origin in the 
mean-field treatment of the local constraints on the pseu- 
doparticle occupation number which is characteristic for 
those theories. 

We analyze the spatial distribution of the impurity- 
induced moment by means of the local Green's function 
G{iuJn;R,R) = -{TrfR{T)fl{0))^„ which is givcu by 



G{iuj„;R, R) ^ G°(iw„) + (5G(iw„; R, R). 



(21) 



Depending upon whether R lies on the A or B sublattice 
the impurity contribution 5G{iujn] R, R) is 

SGiiLUn,R,R) = \A{iLUn;R)\yG°{iLUn), R € A, 

5G{iuJn,R,R) = -\B{iuj,-,;R)\''/G\iuJn), ReB (22) 
with 

A{iujn\ G A) = ^ gl[i^^n] k) exp{ikR), 



X exp(— + ikR). 
Evaluating Eqs. ( p^ ) one finds 

src Tim ^^a{R}_ ( V2R\u;n\\ 

0G[tUJn; R, Rj = ^3:; — An 



(23) 



ttD^ 



D 



X iujn In ^ 7-^ for R £ A, 



SG{iujn;R, R) 



4$B(i?; 

ttD^ 



' D 



X iuJn In"^ for ReB, (24) 
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where the Ki are MacDonald's functions^ and the an- 
gular dependence is determined by 



^ A) = (1 + cos7ri?^) /2, 
$B(-Re B) = (l - {2Rl/R^ - l)cos7ri?:r) /2. 



(25) 



From Eqs. ( |24| ) we determine the impurity contribution 
to the local density of states, dp{ijj\ R) = —\l'KlTtv5G{Lo+ 
iO~^;R,R). It is found that the local moment resides 
predominantly on the B sublattice; for R\uj\ <C J one 
obtains 



dp{uj; Re B) 



2^b{R) 



Siuj) 



1 

Iii{D/\lu\) 2H 
1 



{n/2r + ln' {D/\cj\) 



(26) 



while Sp{u!; i? e A) is negligible. 

Our particular interest is the local spin polarization in 
the vicinity of the impurity as this quantity is needed 
below to derive an expression for the NMR Knight shift 
which we will compare to experimental data. In the pres- 
ence of a uniform magnetic field the spin polarization on 
site R can be obtained from the local spin susceptibility 



d 



(27) 



Evaluating Eq. ([2^) for R\uj\ < J with R e B one ob- 
tains 

6x{T; ReB) = ^ ^^.fiT, R) In'' f (28) 



with a correction factor 



f{T,R)^ixK,{x))\ X 



aV2TR 
ttD ' 



(29) 



which approaches unity for TR <C J. Due to the R~^ 
decay of 5x{T^ R) the spin polarizability is found to be 
strongest on the nearest-neighbor sites of the vacancy. 
Contributions from sublattice A are again negligible. The 
logarithmic correction in Eq. ( |2^ ) is due to the marginal 
character of the bound state: its spectral weight at a 
given lattice site vanishes logarithmically as w ^ 0. In- 
tegrating Eq. ( ^8| ) over R one again obtains the contri- 
bution to the uniform susceptibility 6x{T) = 1/(42^) as 
required. 

Wc now turn to the spinon bond parameter i^R^R+s = 
(ftt+sfR) defined for R e A, which we express in 
terms of the Green's fmiction G{iujn', R, R + S) = 
-(TV/K(r)/],,,(0))<.„: 



Arm+s = -iTj^GiiuJn; R, R + S). 



(30) 



Splitting the Green's function into an unperturbed and 
an impurity-correction term one obtains 



^R,R+S 



5A 



R.R+S 



(31) 



with 



iuJrb 

SAr^r+s = -iTY,5G{ilo^- R,R + d), 



(32) 



where 



G°itLu^;0,d) = B{iuj.^;S), 
dGiitUn, R,R + S)^ A*{iLUn; R) 

X B{iLUn;R + d)/G°{iuj„ 



(33) 



For R = one finds G^{iLUn;0, S) = —5G{iuJn] 0, 5) which 
reconfirms the vanishing spinon bond parameter con- 
necting the impurity site with its nearest neighbors, i.e., 
A0.5 = 0. 



IV. KNIGHT SHIFT 

Experimentally, the spatial distribution of an impurity 
induced moment can be investigated by means of the 
NMR Knight shift K{T) on either the impurity, Cu, or 
in the case of YBa2Gu307 on the Y nuclei. Further infor- 
mation is obtained from the width of the NMR lines and 
the relaxation rates of the nuclear magnetic moment. 

The Hamiltonian describing a nuclear spin I in an ex- 
ternal magnetic field coupled to the surrounding elec- 
tron spins Si is 



Hn ^ InHq ■ I + Am^S^ ■ 7, 



(34) 



where denotes the nuclear gyromagnetic ratio and 
Am is the coupling constant of the isotropic (supertrans- 
ferred) hyperfine interaction. In the presence of a non- 
magnetic impurity the local spin susceptibility that de- 
termines the on-site spin polarization is altered according 
to Eq. (^8|) , yielding for the impurity contribution to the 
Knight shift at distance R from the vacancy 



SK{T; R) 



iAm 



J2Sx{T;R + n), 



(35) 



where the sum is over nearest-neighbor lattice sites of the 
nucleus being probed and j — gfiB with the g factor of 
the Cu^^ spin and the Bohr magneton ps- Taking into 
account only coupling to the most dominant moments 
that reside on the nearest-neighbor sites surrounding the 
vacancy, Eq. (p5h becomes 



5K{T) 



■yAhf 



In 

The correction factor 



— /(r)in-^-. 



(36) 
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200 

T(K) 

FIG. 1. Temperature-dependence of the inverse of the im- 
purity-induced Knight shift 1/5K{T) on the impurity nucleus. 
The solid line is a fit employing Eq. (^6|) to experimental data 
by Ishida et al. for 3% ^^Al-doped Lai.gsSro.isCuOd indicated 
by circles. 



f{T)^{yKM)\ V 



4%/2T 

TlD ' 



(37) 



reduces to unity in the hmit T <^ J. The variable n 
in Eq. ( p6| ) denotes the number of sites that are nearest 
neighbors to the vacancy as well as to the nucleus being 
probed. 

The Knight shift on the impurity nucleus is influenced 
predominantly by the local moments residing on the four 
nearest- neighbor sites of the impurity as described by Eq. 
( ^ ) with n — 4:. Figure ^ contains the calculated T de- 
pendence of 1/SK{T). We set J = 1500 K and g = 2; by 
varying the supertransferred hyperfine coupling constant 
Am the curve is fitted to experimental data obtained by 
Ishida et al.U for the Knight shift of ^^Al performed on 
3% ^'''Al-doped Lai.85Sro.i5Cu04. From the fit we find 
for the supertransferred hyperfine coupling constant a 
numerical value of ^^Aht/Tgcj-ip^gige N ~ ^^■'^ kGs. This 
compares well with ^^^hf/7Ar = 16 kGs found by Ishida 
et al. 

Upon doping the system with impurities, the NMR 
signal of Cu and Y is observed to split into a main line 
accompanied by satellite line(s) as is discussed in the fol- 
lowing. Within the approximation employed in Eq. (|3^ ) 
we assume the impurity-induced moments to reside pre- 
dominantly on the nearest-neighbor sites of the vacancy. 
Resonance signals originating from nuclei that do not lie 
in the vicinity of the impurity are affected only little by 
these moments and contribute to the NMR main line. 
Resonance signals from Cu nuclei on nearest-neighbor 
sites of the impurity are shifted with respect to the main 
line as described by Eq. ( p6| ) with n — 1 and constitute a 
single Cu satellite line. Resonance signals from Y nuclei 
located on nearest-neighbor and next-nearest neighbor 



sites of the impurity are seen as two satellite lines cor- 
responding to n = 2 and n = 1, respectively. These Y 
satellite lines have been observed experimentally by Ma- 
hajan et alB in 1% Zn-doped YBa2Cu3 6.64- 

We now discuss briefly the NMR line-broadening which 
results from the inhomogeneous distribution of local 
Knight shifts in the presence of many vacancies. Accord- 
ing to Eq. ( p8| ) moments induced by a vacancy located 
on the A sublattice are found on all sites of sublattice 
B with magnitude oc R~^, R being the distance from 
the vacancy. As a result the nuclear resonance signals 
contributing to the main and satellite lines vary with 
distance from the impurities which can be observed as 
a broadening of the line. At a finite concentration c of 
vacancies the linewidth AH can be estimated from Eq. 
( ^ ) , assuming an average distance between vacancies of 



AH ^ "fAM c 
Hq jn 27rr 



In" 



D 
f' 



(38) 



V. IMPURITY INTERACTION 

In this section we study the interaction between two 
vacancies, first by considering two isolated impurities. 
The case of a finite concentration of impurities is treated 
subsequently. 

The spinon Green's function in the presence of two 
impurities located at sites and R is 



(39) 



Employing scattering theory the two-impurity correction 
5Git{iujn) can be expressed as 

SGH,{iUn) = 

1 ^ / G"(*c.„;r,0) \ f G«(zc.„;0,r) 

n2^[ G^iiio^- r,R))^ 1^ G%^u;„■, r) 

(40) 

with the unperturbed Green's function G'°(iw„; R, R') = 
— (Ir/ji(T)/j^/(0))°^ and the scattering matrix 



T21 T22 



(41) 



with 



Til =r22 = -l/G(zw„,fi,i?), 

T12 = G°{iLu„; 0, R)/ [G{iu;„; R, R)G°iiu;„)] , 

T21 = G"(zc^„;H,0)/ [GiiUn;R,R)G'>{iUn)] . 



(42) 



The diagonal and off-diagonal elements of T [i{iujn) de- 
scribe successive scattering by one impurity and by the 
two impurities, respectively. From Eqs. (pO|)-(^) one ob- 
tains 
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5GR{iuin) 



diuj, 



■ln(G°(iw„)G(zc^„;H,i?)), (43) 



where G{iu)n] R, R) is given by Eq. (pi]). 

By means of Eq. (Ba) we evaluate the two-impurity 
contribution to the density of states Spji{u!) = 
— l/TTlm6Gii{uj + 0"''). For R ^ B one finds a sphtting 
of the resonance energy level 



SpResiio) ^6{iu^ J{R)) + d{uj + J{R)) 



(44) 



where the parameter J(-R) controlling the level splitting 
is given by 



J{R) 



D^2<^b{R) _i 
2R ^ 



ttR 



V^'^BiR) 



(45) 



The level splitting results in an attraction between the 
two impurities due to the formation of a singlet of the 
induced moments. For the shortest distance between va- 
cancies, R = 1, Eq. ( |45| ) gives J{R) « J as expected: 
J is the only energy scale of the model under consider- 
ation. In the present mean-field theory no attraction is 
found between impurities lying on different sublattices, 
i.e., R G A. These findings are also consistent with ex- 
act diagonalization data of Bulut et alS3 who found a 
binding energy of —0.56 J for nearest-neighbor static va- 
cancies on a 4 X 4 lattice, while the interaction between 
vacancies separated by R is negligible if Rx + Ry is an 
even number. We interprete J(R) given by Eq. ( ^5|) as 
being the antiferromagnetic exchange coupling of impu- 
rity moments. Indeed, the impurity contribution to the 
spin susceptibility calculated from Eq. (H^) is 



6x{T » J{R)) 



1 

2r' 



5x(r« J(i?)) = -exp 



J{R) \ 



(46) 



which is precisely that of two 1/2 spins coupled by an 
interaction Hiat = J{R)Si ■ 82- 

From Eq. (^) the interaction potential between two 
moments is found to fall off slowly as R~^ with distance 
between the vacancies. In a real system with finite im- 
purity concentration c, however, one expects this inter- 
action to be screened at large distance by the presence 
of other impurities. To account for this effect we intro- 
duce a finite self-energy S(a;) determined by means of a 
coherent-potential approximation (CPA): 



Y.{uj)G°{lu -1:{lo)) +c = 0. 



(47) 



Neglecting the frequency dependence of we approx- 
imate the self-energy by S(cj) « I](w ^ 0) = — zISqI, 
where K fulfills 



D 



c = 0. 



(48) 



Reevaluating the exchange coupling parameter J{R) 
one finds an exponential cutoff in the two-impurity inter- 
action potential 



J{R) (X exp {-R/R^^) 



(49) 



where i?av = -^/l/c is the average distance between im- 
purities. For small distances i? <C J one recovers again 
Eq. (^). The present CPA treatment gives a qualita- 
tive description of the screening of the interaction po- 
tential of two distant impurities, but it neglects impor- 
tant effects such as the formation of spin singlets among 
closely spaced moments; a complete analysis of the many- 
impurity problem should take this effect into account. 

Having the interaction scale between moments we can 
give an estimation of the relaxation rate 1/Ti of nu- 
clear magnetization resulting from fluctuations of the 
impurity-induced moments in the vicinity of the vacancy. 
To be specific let us consider the nuclear-spin relaxation 
rate 1/Ti at the impurity site (obtained, for instance, 
from the 2^A1-NMR signal in an experiment by Ishida et 
alB) With the supertransferred hyperfine interaction of 
Eq. ( ^ ) the nuclear relaxation rate is 

1 



2 .2 



M 



(50) 



SS' 



where r is the correlation time of the local moment and 
h = I. The amplitude of spin fluctuations on nearest- 
neighbor sites of the impurity is given by the equal-time 
spin-correlation function 



SS' 



[Ss ■ Ss') 

-ST^ 'Y6G{iu„ + iiym;S,S')dGiiuj.^;S',S) 



5,, (51) 

where SG{iujn',S,^') = —B*{iuJn;6)B{iuJn;5')/G^{iuJn), 
Vm — ^mnT are bosonic frequencies, and z = 4 is the 
number of nearest neighbors. We assume the correlation 
time T in Eq. ( |50| ) to be dominated by the exchange in- 
teraction among moments. In principle, r has a certain 
distribution due to the random spreading of impurities. 
We do not go into this delicate issue and roughly approx- 
imate T « 1/ J(-\/2_Rav)] , where J(-\/2i?av) is given by 
Eq. (H) and V2i?av = is the average distance be- 

tween impurities located on different sublattices. One 
then obtains 



V2 



In- 



(52) 



Assuming an impurity concentration c = 3%, a temper- 
ature T — 100 K, and using the value of the supertrans- 
ferred hyperflne coupling constant ^hf found in the pre- 
vious section, Eqs. (|^)-(|5|) give 1/Ti w 0.29 ms^^. We 
compare this, value to l/Ti — 0.43 ms~^ observed by 
Ishida et al.U for 3% ^''Al-doped Lai.85Sro.i5Cu04 and 
flnd it to be of the same order of magnitude. 
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VI. KONDO SCREENING OF LOCAL MOMENTS 

The two-dimensional Heisenberg model has a Neel- 
ordered ground state and it was implicitly assumed above 
that the spin-gap disordered state is stabilized by mobile 
holes. Naturally, the question then arises how a local 
moment induced by a static vacancy will be affected by 
these charge carriers. Once the metallic state with fi- 
nite Fermi surface is formed one usually expects Kondo 
screening of the local moment at low temperature. As 
an exception, a spinon bound state on a vacancy can 
coexist with metallic conductivity if spin-charge separa- 
tion occurs with a gap in the spin-excitation spectrum. 
Underdoped cuprates, however, seem to have ungapped 
pieces of their original (large) Fermi surfaces even in a 
"spin-gap" phaseE2l. Therefore the bound spinon on the 
vacancy is screened by gapless Fermi-surface excitations, 
although the Kondo temperature is expected to be low 
due to the reduced density of states. The characteris- 
tic Kondo-energy scale can be roughly estimated in the 
following way. 

We introduce fermionic field operators ipia = Sa^{Ri) 
which annihilate a spinon Sa- in a bound state with wave 
function ip{Ri), Ri being the distance from the spinless 
vacancy. Due to the hopping term in the t- J model there 
is a finite amplitude tip{Ri) of mixing of the bound spinon 
with extended states: 



H = t 



i,S,(T 



H.c, 



(53) 



The bosonic operators h] and c|^^ ^ in Eq. (^) create 
a "holon" and a conduction electron, respectively. In a 
slave-boson formulation, the latter is a product of the 
lattice spinon and holon operators. 

At finite doping, x — ^ 0, the Hamiltonian (|52 

allows for a localized spinon to escape the impurity which 
leads to a broadening of the bound state. The width of 
this resonance controls the low-temperature behavior of 
the local spin and is given by the Sa fermion self-energy 

^%u;)=t'Tj2 E ^*iR^MRJ)Gt{R^,) 

V ij.SS' 

xGl,,{R,+s-R,+s'). (54) 

At low temperatures the holons are almost condensed, 
and the holon Green's function can be approximated by 
« {x/T)5^fi. Furthermore, we keep only the nonoscil- 
latory on-site terms that present the leading contribu- 
tions to Eq. ( p4[ ) and employ the normalization condition 
X^Jv'iP = 1- Then the width of the Kondo resonance is 



lml]" (0) w TTzxt^N^iO), 



(55) 



where N'^{0) — I/tt Im 6*^(0-1 — 0) is the on-site den- 
sity of electronic states on the Fermi level. 

In the overdoped t-J model with a large gapless Fermi 
surface the density of states is controlled by the inverse of 



the full bandwidth, A^^(O) « l/{2zt). In the underdoped 
pseudogap regime, xt <C J^ the density of states A^'^(O) 
is strongly reduced which is due to the fact that the un- 
gapped areas of the Fermi surface contributing to iV'^(O) 
are only small. We simply assume A^'^(O) w xj{2zt) and 
estimate the Kondo temperature using Eq. (|55| ) as 



Tk{x) 



-xH. 



(56) 



For a typical value of x w 0.1 and t w 0.4 eV this 
gives Tk ~ 70 K. Tk{x) quickly increases with x, caus- 
ing the vacancy-induced moment to disappear in the opti- 
mal and overdoped regime. This magnetic-nonmagnetic 
Kondo crossover upoa hole doping was recently employed 
by Nagaosa and LeeEl to explain the unusual doping de- 
pendence of the residual resistivity of cuprate supercon- 
ductors. 

One may wonder why the exchange parameter J does 
not enter in Eq. (|5^). In fact, the main role of the Heisen- 
berg term in this quite unusual Kondo-like behavior of 
a spinless vacancy is to produce a spin pseudogap hence 
creating a local moment. It is the hopping term that is 
then responsible for converting the spinon bound state 
into a resonance of finite width. 



VII. CONCLUSION 

In summary, we have studied the localized magnetic 
states induced by static spin vacancies in underdoped 
high-Tc cuprates. Starting from the two-dimensional 
Heisenberg antiferromagnet we employ a drone-fermion 
mean-field theory describing a spin liquid in the spin- 
gap phase. The impurity induces a spinon bound state 
residing predominantly on the nearest-neighbor sites of 
the impurity. We have calculated the impurity-induced 
Knight shift and found it in mod agreement with ex- 
perimental data of Ishida et alu for Al-doped cuprates. 
The calculations do also account for the appearance of 
satellite peaks in NMR measurements on Y observed by 
Mahajan et alu in Zn-doped cuprates. Two impurities 
are found to attract each other due to the formation of a 
singlet of the induced moments. In the presence of a fi- 
nite concentration of impurities the interaction between 
moments is screened at distances larger then the aver- 
age separation between impurities. We estimate the con- 
tribution of local moments to the NMR relaxation rate 
which is found to be in reasonable agreement with exper- 
iment. Finally, an estimation of the Kondo temperature 
of screening of the vacancy-induced moments by charge 
carriers is given, below which a nonmagnetic impurity 
behaves as a sharp resonance at the Fermi level. 
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